Distinguishing Spin Relaxation Mechanisms in Organic Semiconductors 
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A theory is introduced for spin relaxation and spin diffusion of hopping carriers in a disordered system. For 
disorder described by a distribution of waiting times between hops (e.g. from multiple traps, site-energy disorder 
and/or positional disorder) the dominant spin relaxation mechanisms in organic semiconductors (hyperfine, 
hopping-induced spin-orbit, and intra-site spin relaxation) each produce different characteristic spin relaxation 
and spin diffusion dependences on temperature. The resulting unique experimental signatures predicted by the 
theory for each mechanism in organic semiconductors provide a prescription for determining the dominant spin 
relaxation mechanism. 



Spintronics in organic semiconductors fTl provides dra- 
matically different regimes than are common for inorganic 
semiconductorslJl O or metalsH |5l, due to the very low 
mobilities and conductivities of organic materials, as well 
as the ubiquity of constituents with low atomic number (and 
thus often weak spin-orbit coupling (SOC) and long spin re- 
laxation times). A major focus in organic spintronics has 
been on organic spin valves |6 -10| , within which spin po- 
larized carriers are injected from one magnetic source con- 
tact into a non-magnetic organic spacer layer which they tra- 
verse before being collected by a magnetic drain contact. If 
the spins maintain their polarization over the length of the 
spacer region, then parallel and anti -parallel electrode magne- 
tizations produce different resistances and lead to spin-valve 
phenomenology. The spin-valve effect requires the preser- 
vation of spin within the spacer layer; therefore the reported 
long spin lifetimes in organic semiconductors increase organic 
spin valves appeal. Any spin-based application beyond a spin- 
valve should also require an understanding of the rate of decay 
of the spin polarization in organic semiconductors. The domi- 
nant spin relaxation mechanism is controversial llTl fTTHTSl but 
appears to be driven by either SOC or the hyperfine interac- 
tion (HFI) during hopping, or associated with an intra-site re- 
laxation (ISR) process that may include spin interactions with 
phonon modes or the multiple nuclear fields at a molecular 
site (independent of hopping). Only a few theoretical investi- 
gations lfT2lfT4l have examined spin relaxation including dis- 
order, even though the interplay between spin and charge dy- 
namics is well-known for inorganic semiconductors! 16 1. 

Here we provide a unified theoretical description of spin 
lifetimes in organic systems from SOC, HFI or ISR, based on 
a continuous -time random-walk theory ifTTlfTSl . that predicts 
dramatically different spin lifetime dependencies on temper- 
ature as well as previously unrecognized analytic dependen- 
cies on other physical properties of the organic semiconduc- 
tor; these results are summarized in the upper half of Table 
I. We introduce disorder in two scenarios: (1) by allowing 
hopping caiTiers to be captured by traps at a single energy 
or (2) at a wide aiTay of energies. The second scenario uses 
a hopping-time distribution that lacks a first moment (a so- 
called heavy-tailed wait-time distribution (WTD)) and also 
describes charge transport when energetic disorder of a cer- 
tain type is present that produces time-dependent mobility 



(so-called 'dispersive' transport). These heavy-tailed WTDs, 
which at long times scale as f^'^" where a is a character- 
istic parameter of the WTD, provide the phenomenology as- 
sociated with each spin relaxation mechanism (SOC, HFI, or 
ISR) during dispersive transport. Dramatic qualitative differ- 
ences emerge between HFI and SOC spin relaxation times (Xj) 
and spin diffusion lengths (£s), which positions our theory to 
predict experimental observations to differentiate between the 
two; e.g. our theory provides supporting evidence that SOC 
dominates in Alqs whereas HFI dominates in MEH-PPV. In 
addition, HFI, SOC, and ISR are incorporated simultaneously 
within a single theory, which illuminates cross-over behavior 
from one mechanism to another Our results for low disor- 
der, reported in the non-dispersive section of Table I, agree 
with previously-obtained analytic results for SOC lfT4l and 
HFI ||T9 l spin dynamics. We further agree with Ref.TTs result 
for the spin diffusion length in the fast hopping HFI dispersive 
regime for the special case of a = 1/2. By varying tempera- 
ture, HFI strength, and SOC, the entries of Table I change in a 
unique manner; our theory will aid experimental efforts to dis- 
tinguish which relaxation route occurs. As our theory uses an 
arbitrary WTD, it is applicable to a wide variety of other types 
of disorder and materials, such as spin transport in amorphous 
silicon or colloidal quantum dot films. 
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TABLE I. Spin relaxation rate (x^') and spin diffusion length (£,) 
dependencies on width of the Gaussian HFI distribution (a), SOC 
(y), hopping rate (ko), and ISR (F). ko, T, and a depend on temper- 
ature, and thus allow access to different regimes in the same sample. 
Deuteration affects a and elemental composition affects j. Fast hop- 
ping HFI considers only short-timescale spin relaxation/diffusion. 
Spin decay in the dispersive regime is algebraic in time. 

Theory with general WTD - The dynamics of a classical 
spin S in an arbitrary precessional field CO, influenced by an 



ISR process with a phenomenological rate F, is 

dS{t) 



dt 



:mx5(f)-r5(f) ==(£i-ri) •5(f) = ((o£i-ri) -5(0, 

(1) 

where ii is a skew-symmetric matrix |20|, and is solved by 
5(f) = e-rfgflr . ^^ = g-rr^j-j-) . ^^^ where 5o is the initial spin 
vector and 

R(f) = i+sincofi2 + 2sin2— Q-Q. (2) 

Here we consider fields from two sources: a hyperfine field 
(fi)/,/) and a spin-orbit field (y). Spins at different sites which 
begin in the same polarization state experience different hy- 
perfine fields, which causes partial dephasing of the spin en- 
semble fJTl; polarization is further lost when the spins inco- 
herently hop and experience different hyperfine fields at other 
sites lf22l . Hopping also affects spin polarization through the 
spin-orbit interaction; spin-flips are possible as the carriers are 
not in pure spin states [14 , 23 , 24 1 . Figure [TT a) schematically 
indicates the fields a spin might experience over a time in- 
terval due to hyperfine and spin-orbit fields. Lastly a spin at 
any site may interact with phonon modes which leads to ISR 

To describe the evolution of an ensemble, we simplify the 
single-spin dynamics by assuming: 1) hyperfine and spin- 
orbit rotations are independent and isotropic in space, 2) hy- 
perfine rotations occur discontinuously at each hop 1261 (as 
represented in Figure flTa)), 3) spin-orbit rotations of magni- 
tude Y (also treated phenomenologically) are independent of 
the duration of time spent at a site (they only occur at hops) 
and 4) ISR is independent of hopping. Hyperfine fields are 
distributed according to a Gaussian distribution with width a - 
a phenomenological parameter assumed to be on the order of 
1 mT. The random direction of spin-orbit fields results from 
the random orientation of localizing molecular sites lfT4l . 

After averaging Eq. l2]over the hyperfine field distribution. 
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whose elements are found in II20II . The spatially averaged 
spin-orbit rotation matrix is Ry = 1 — (4/3) sin^(Y/2) 1. Typi- 
cally, Y is assumed to be a small angle such that Ry ss 1 — ^ Y^ 1 ■ 
We now consider the spin ensemble to be hopping from 
site to site where the WTD is arbitrary; the time evolution of 
the spin polarization can be calculated and the spin relaxation 
time can be extracted - numerically in general, but analyti- 
cally in certain cases to be discussed. Since the HFI oper- 
ates between hops and the SOC affects spin at the hopping 
event, their rotation matrices enter into the polarization func- 
tion in different ways which preclude interference effects. In 
this sense, HFI is analogous to the D'yakonov-PereF 1271 and 
SOC is analagous to the Elliott- Yafet [28) |29l spin relaxation 
mechanisms for inorganic semiconductors. 



The rotation accrued during the wait-time is Rj, which is 
'Rhf when HFI is present and 1 when it is absent. An angular 
rotation R/, = Ry occurs from the hopping event when SOC 
is present but is 1 when SOC is absent. The third spin loss 
avenue, intra-site spin relaxation, simply appends exp (— Ff ) 
to Rj as seen in Eq. ([iji. With these different factors in mind, 
we write the R due to all single hops in the spin ensemble as 

Ri(f)=R,,yJ^'R,(f-f')<I>(f-f')e"^<'"''*(R.(f')¥(f')e"'~'')^f'- 

(4) 
The hopping rotation matrix is factored out of the integral 
since it is independent of time. x|/(f ) is the WTD from which 
the wait-time f is drawn. It appears once since one hop occurs. 
<I>(f) is the survival probability: the probability that the wait- 
time at a site exceeds f: 4>(f) = /,°°\|/(f')(if' [30|. The quanti- 
ties Ro(f) = Rs(f)4>(f)e-'^' and R'o(f) = R,,(f)x|/(f)e"'^' permit 
the compact expression Ri (f ) = R/, /q Ro(f - f')Ro(f')*'- The 
procedure can be continued indefinitely for an arbitrary num- 
ber / of hops yielding the recursive equation 



R,+i(f)=R/, / R[)(f-f')R/(f')^f', 
Jo 



(5) 



which has the form of a convolution; R/(i) — RJjRo(i + 

F)Rq [s + F) in Laplace space. We build the polarization func- 
tion from the different number of random rotations ll22l [3T1 : 



P(f)=P(f)-So = £R,(f)-So, 



(6) 



where R/(f) So is the polarization after I jumps occurring 
between time and f. The polarization can be calculated in 
Laplace space by summing the geometric series Eq. l6]to ob- 
tain our main equation: 



P(i)=Ro(.9 + F)[l-RAR;(i + F)]-i-S, 



(7) 



We have left the polarization as a vector since when an applied 
field is present, the rotation matrices need not be isotropic and 
the polarization in general decays with different longitudinal 
and transverse spin relaxation rates. 

An exponential WTD, \(/(f) = fcexp(— ^f) characterized by 
an average hopping rate k, provides a clarifying example of 
the theory, as RQ(f) = feRo(f) and the Laplace transform of 

Ro(f) is especiaUy simple: jQRo{t)e^^'dt — Rs{s + r + k). 
Eq. ITlbecomes 



?(,?) == R, (i + F + fc) [1 - fcR,,R, {s + T + k)]-'-S 



(8) 



With HFI present, the polarization cannot be solved analyt- 
ically in the time domain except in special cases. This po- 
larization function has been studied in the context of muon 
spin rotationQH |3ll ES). For only SOC (and F == 0), the 
polarization in the time domain is calculated exactly to be 
P{t) = e^'^l^'", which is in qualitative agreement with Yu's 
recent determination U4, ^ . If Y = 0, then P{s) ^l/{s + r) 



(i.e. P{t) — exp(— r?)) can be shown from the relation be- 
tween <l> and \j/. Ref. 20 contains a general expression for Eq. 
(|7]i in the absence of hyperfine fields. 

Traps: simple models - Before examining systems with re- 
alistic WTDs, three instructional models are investigated that 
have analytic solutions. The easiest model to study is that of 
a single trapping level. There are two ways to obtain nearly 
equivalent results: 1) incorporate explicit traps by forming a 
modified renewal equation in the spirit of Eq. dTl and 2) use 
a bi-exponential WTD to describe normal hopping and trap 
release. The first approach was used in a limited manner by 
Kehr and Honig 12211 . We have generalized their approach 
to include SOC but only report the results here in Figure [T] 
(b) (black lines) where analytic solutions for the polarization 
function are achievable if a = 0. 
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FIG. 1. (Color online) (a) Schematic of the time-dependent fields a 
single spin experiences from HFI and a constant (y) SOC. For (b), 
(c), and (d): k, the average hopping rate, is 100, a = and F = so 
these results are for SOC only, (b) Polarization as a function of time 
for explicit trap (black) and bi-exponential (red dotted) models. The 
rate of leaving a trap, kQ, for (a), ko = 10 and (b), A;o = 0.1 where 
each pair is labeled by k\, the rate to be captured by a trap. Spin 
relaxation decreases as k\ increases (i.e. as trap number increases). 
The two models are connected through the relations b = k\l{k^k\) 
and (3 = ki^jk. Multi-exponential model in (c) with varying SOC 
labels and in (d) with different b. Slopes of long time power law 
dependence are indicated. 

The second approach accounts for traps by using the bi- 
exponential WTD: \(/(f) == (1 +Z7)-'A:e-^' +/7p(l +/7)-iA:e-P'^' 
with < {^,P} < 1. We consider b to be controlling the 
importance of traps. When h is small, traps become incon- 
sequential as seen from the fact that the single exponential 
WTD emerges. The dimensionless quantity P determines the 
reduction of the trap release rate compared to the normal hop- 
ping rate, k. It is important to mention that spin relaxation re- 
sults within the two-state trap model and bi-exponential WTD 
are always fit well with exponential decay times (aside from 
short time dynamics). Figure[TJb) shows an example of SOC- 
only spin relaxation (red-dotted lines) and compares the ex- 
plicit and bi-exponential trapping models. Since spin decay 
is produced by hops, long release times (l/^o) promote long 



lived polarizations as shown in Figure [Tfb). Likewise more 
traps (larger k\ ) reduces overall number of hops and spins live 
longer 

Heavy-tailed WTDs have successfully been used to de- 
scribe charge transport in disordered systems 1,17,1 . The suc- 
cess of the bi-exponential model suggests an easy way to ex- 
tend the model to a diverse array of trap frequencies. We 
construct a multi-exponential WTD from the bi-exponential 
model as follows: we allow for more trap release times 
by taking p — > p' with j being a positive integer. Trap- 
ping probabilities are also generalized \a h ^ V . The re- 
sulting WTD has been previously studied 0414361 : \|/(f) — 
^"'(1 -^)I7=o^^+'P^'te-P''='. For the condition /; < p, this 
WTD does not have an average wait-time which distinguishes 
it from the bi-exponential WTD. At long times, the algebraic 
form is found to be \|/(f) -- f"""' where a = In/j/lnp M30II34I . 
Figure [TJc, d) show different scenarios with this heavy-tailed 
WTD. Axes are log-log to emphasize the algebraic decay as 
opposed to the exponential decay in Figurefltb). 

Mutiple-trapping - We now examine a realistic WTD — 
one due to trapping levels distributed exponentially and for 
HFI and ISR as well as SOC. It is now thought that hop- 
ping within an energetic distribution of hopping sites is much 
like what occurs within the trapping model that is consid- 
ered here ||37l 1381 ; essentially in each case the deep en- 
ergy states lead to dispersive transport Ii391 . In general 
the mutiple-trapping WTD is written as IJOJ [41] \|/(0 = 
/_''„£/eg(e)^(e)exp(— ^(e)f), where g(e) is the density of 
states and k{e) ~ koe^'^'^^ . Hartenstein et al. showed fSSl 
that this multiple-trapping WTD agrees well with the Monte- 
Carlo-simulated hopping WTD. 

At present we consider only an exponential density of 
states, and note that the end results are also valid for a wide 
Gaussian density of states (a must be changed to 1/2); it is 
advantageous to write the WTD in Laplace space: 
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(9) 



where x — e/ksTo. The long time behavior of the distribu- 



tion is algebraic \\i{t) ■ 



,-i-a 



where a=T/To\40\ (here T is 



the temperature and To is a parameter characterizing the band- 
width of the states involved in transport). For SOC, the inte- 
gral can be computed so the inversion techniques of Ref. 42] 
can be utilized to solve Eq. (JTl). However, the form of the spin 
relaxation for both SOC and HFI can be deduced from the fol- 
lowing physical arguments. From continuous-time random- 
walk theory |30|, the number of hops scales as « ^ (^of)"- 
Unsurprisingly for SOC, the polarization falls with increasing 
number of hops; it can only do so given spin-flips which hap- 
pen with probabihty ^ y^. So in general P{t) ^ {/{koff-^"-)"- 
for a heavy-tailed WTD with SOC. The effect of HFI is more 
subtle. When hopping is slow {ko < a), even free hops lead 
to rapid spin randomization - on the order of the hopping time 
and independent of the strength of the HFI. When hopping 
is fast and traps ineffectual, normal hops (a > 1 in Figure l2]) 



lead to very little polarization loss since motional narrowing 
occurs. When deep traps are present (a < 1) then fast hop- 
ping entails that deep traps are quickly populated; the result- 
ing polarization behaves similarly to the Kubo-Toyabe polar- 
ization function |21 1. This behavior is characterized by two 
spin relaxation times: at short times, XsO ^ 1 (blue diamonds 
in Figure [2]), and at long times a much longer relaxation time 
exists (red triangles) that approaches zero as the trap density 
increases. The fast relaxation at short times is strikingly dif- 
ferent than free hopping in that the rate is independent of k(). 

This reasoning can be proved for SOC - the long time polar- 
ization is determined analytically to be algebraically decaying 
as P{t) ^ t 
1/e) is 



201 . The spin relaxation rate (rate to decay to 



1 f /tof /"e-'/"ri/« a<l 
ti [ kof- a > 1 



(10) 



and is also displayed in Figure l2] (solid black line) where the 
role of disorder is apparent. / has further a dependence but 
is weak when a not near one. Since 7b is a measure of the 
disorder, larger disorder increases T^ which agrees with our 
single trap analysis. Alternatively, this analytic result can 
be checked for the exponential density of states and another 
heavy-tailed WTD (multi-exponential) by finding where the 
numerical polarization functions decay to 1/e (with a de- 
fined appropriately). The HFI within the multiple-trapping 
picture poses a much more difficult problem analytically and 
numerically. However it is quite accessible when the multi- 
exponential WTD is utilized. By analogy with SOC, those 
results are carried over to the multiple-trapping picture with 
the exponential density of states. 

For dispersive transport, a generalized diffusion coefficient 
is defined as Ka = r^A:"/6r(l + a) where r is the root-mean- 
square hopping distance 1301 . This leads to 4 — \/6Kaif- 
In the non-dispersive regime, the diffusion coefficient is D = 
r^feo/6 which entails is — V^DZs- As a confirmation of our 
method, we remark that l^ (when using a short relaxation 
timescale) for fast-hopping-HFI agrees with the analysis of 
Ref 12 for the special case of a = 1/2. One remark- 
able feature is that SOC's 4 has the same form whether the 
transport is dispersive or non-dispersive. Whereas the SOC 
spin relaxation rates decrease due to multiple trapping, mo- 
bility decreases as well so it is not obvious whether 4 in- 
creases or decreases in the dispersive compared to the non- 
dispersive regime. We find that the 4 does not change sig- 
nificantly across the two regimes so that Ref jJJ and 1331 s 
Efros-Shklovskii hopping prescription for Alqs ||9| remain 
valid. Slow hopping HFI relaxation does not explain the data 
since the measured £s is always greater than the typical hop- 
ping length. Alqs is somewhat unique in that its SOC is 
large (y « 0.03 |[33l ) compared to other frequently-studied 
organic semiconductors. For instance SOC in MEH-PPV is 
100 times weaker than in Alq3 ||33]| . This explains why HFI 
has a large effect on spin preservation in the polymer DOO- 
PPV Iil31 . When deuterated, the spin relaxation rate change 



is nearly described by the relation found here: zj^ °^ a where 
fl sa 0.1 — 0.6 rad/nsmT |20|. The change in 4 also agrees if 
a is taken to be near unity; changes in 4 are more difficult to 
gauge since it is not a strong function of a. 

We conclude, by examination of the elements of Table I, 
that the spin relaxation mechanism in assorted organic semi- 
conductors can be probed by altering the HFI (through deuter- 
ation), the hopping rate (through temperature), and SOC 
(through molecule choice). This theory should provide a 
framework for future experimentalists to use for determining 
spin processes in organic semiconductors. In addition, the the- 
ory should apply to other disordered systems in which trans- 
port is dispersive, such as amorphous inorganic semiconduc- 
tors and colloidal quantum dot films [43,1 . 
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FIG. 2. Numerical SOC spin relaxation rate (black solid symbols) 
and HFI spin relaxation rate (red triangular and blue diamond sym- 
bols) as a function of a. SOC calculations use multiple trapping 
WTD with exponential DOS while HFI uses the multi-exponential 
WTD. Solid black line is fit using Eq. ifTol with one multiplicative fit 
parameter. Blue diamonds depict the short time spin relaxation. Red 
triangles represent the longer time slower spin relaxation. Vertical 
gray line separates regions of algebraic and exponential spin decay. 
ko/a = 1000, p = 0.1, and F = 0. 
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